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Abstract 

We study the dynamic fluctuations of the soft-spin version of the 
Edwards- Anderson model in the critical region for T T^. First we 
solve the infinite-range limit of the model using the random matrix 
method. We define the static and dynamic 2-point and 4-point correla- 
tion functions at the order 0{1/N) and we verify that the static limit 
obtained from the dynamic expressions is correct. In a second part we 
use the functional integral formalism to define an eflfective short-range 
Lagrangian L for the fields SQf^ {ti,t2) up to the cubic order in the 



series expansion around the dynamic Mean-Field value (5"^(ii,t2)- We 
find the more general expression for the time depending non-local fluc- 
tuations, the propagators [{SQ'^^ {ti,t2)SQ'^^ {t-i,t4))^]j, in the quadratic 
approximation. Finally we compare the long-range limit of the corre- 
lations, derived in this formalism, with the correlations of the infinite- 
range model studied with the previous approach (random matrices). 
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1 Introduction 



The static and dynamic Mean Field (MF) theory of Spin Glasses (SG) sys- 
tems for T > Tc is well defined and understood. This theory has been studied 
through different approaches all consistent among themselves. Many impor- 
tant results concerning the equilibrium static properties of SG have been de- 
rived using the replica method [||]. Sompolinsky and Zippelius 0, |^ 
studied a soft spin version of the Edwards- Anderson model with the dynamic 
formalism, avoiding the replica trick. They defined a Langevin dynamics on 
the system and analysed the infinite range limit where the MF solution is ex- 
act. The static limit derived from the dynamic expressions is in agreement with 
the static prevision obtained with replica. Moreover, dynamic characteristics 
of the model have been well defined. 

Unfortunately the behaviour of short-range SG system is not clear yet. 
There are different analytic works and simulations about the Ising and Heisem- 



berg SG in finite dimensions, 0, 0, [jTO|, |TT|, that do not agree with each 
other. 

On the other hand, in order to study the corrections to the MF behaviour of 
the Green functions one can use the Renormalization Group theory and the 
e-expansion. Chen, Lubensky [|T^], [|T7| and Green |]TB[ studied a static Ising 
model in d = 6 — e dimensions (6 is the upper critical dimension) with the 
Replica method and they found the corrections to the second order in e for the 
critical exponents. 

In this work, we want investigate the behaviour of the dynamic fluctuations 
of the short range SG in the critical region for T T+. We study the soft- 
spin version of the Edwards-Anderson model that evolves through Langevin 
dynamics. We adopt the same procedure that Sompolinsky and Zippelius 
used in and we manage to explicit write the time- dependent propagators 
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[{5Q^ {ti,t2)5Qj {tz,t4))^j for any value of ti,t2,^3,^4 in the critical region, 
while in they were defined only in the limit of two complete time separation. 
They are the elementary building blocks in a renormalization group calculation 
expanded in e = 6 — and can be used for future studies of the dynamic 
effects of higher-order terms (the cubic interactions). For evaluating these 
propagator we define an effective Lagrangian of the fields (which 
represent the fiuctuations around the MF order parameter value Q'i^{t,t')) 
through the functional integral formulation of the dynamics. In order to have 
a comparative term order by order in perturbation theory, we solve the infinite- 
range limit and the 0(1/N) corrections of this dynamic model, by using the 
distribution of the eigenvalues of the random interaction matrix. We verify 
that the expressions derived with the two different and independent methods 
are consistent each other. 

The aim of this work is to pursue the study of the quadratic fluctuations 
of the soft-spin model in the general case, without having recourse to the 
Glauber model (hard-spin limit) [0. We deflne the quadratic fluctuations as a 
perturbative series in g (the coupling constant of the fourth order term of the 
soft-spin Lagrangian) which we succeed in resuming and therefore in obtaining 
a g independent expression that could be directly used for further diagram- 
matic expansions of the theory. The results are qualitatively different from 
those obtained by Zippelius and this may be related to the approximations 
done in going to the hard case. 

This paper is organized as follows: in section 2 we deflne the theory; in 
section 3 we flnd the form of the quadratic corrections to the MF limit (order 
0(1/A^), where N is the sites number) using the diagonalization of the inter- 
action random matrix; flnally in section 4 we flnd the general propagator for 
an effective short range Lagrangian of the fleld 5Q°'^{t, t') (fluctuations around 
MF limit); in the conclusion we present the possible development of this work. 
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2 Definition of the model 



Let us consider the soft spin version of the Edwards- Anderson model given by 
the Hamiltonian 

where the sum (ij) is over z nearest-neighbor sites and the couphngs Jjj are 
quenched random variables with the following distribution: 

m.) = ^^^exp(-.4./2). (2) 

A purely relaxational dynamics can be described by the Langevin equation: 

^° dt ~ ds,{t) ^"^^ 
^i{t) is a Gaussian and white noise with zero mean and variance 

i 

which ensures that locally the fluctuation-dissipation theorem holds. 

In the MF theory the physical quantities of interest are the average local 
(at the same point) correlation function 

C{t-t') = [{s,{t)s,{t')),]j, (4) 

and the average local response function 

G{t-t')^^[{s,mi{t'Mj, (5) 

where the angular brackets refer to averages over the noise and square brack- 
ets over quenched disorder J^j. 
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Beyond the MF approximation we evaluate the non local fluctuations that 
are non vanishing 

where a,/3,7,5 can take the values 1 or 2 being cj)} = and 0f = Sj. We 
shall see that these quantities represent the propagators of our theory. We will 
focus on the properties of the small frequency and small momentum of the 
propagators that have a critical slowing down near T^. 

In this formalism the time- depending spin-glass susceptibility is 

XsG{^-j■,tl-t3,U-h) = j^[{si{ti)Cj{t3))i{sj{U)Ut2))i]j = Gl^^\t-j;h-h,U-t2), 

(7) 

where the Gq^'*''^ functions represent the fluctuations in the limit |ti — ^2! ^ 00 
and 1^3 — t^l ^00. 

3 Mean Field limit and quadratic fluctuations 
by diagonalization of the interaction 
matrix 

The theory defined with the Hamiltonian (|l|) where the Jij are quenched ran- 
dom interactions can be solved by using the general properties of random 
matrices. It is convenient to apply this approach to an infinite range model 
{z = N and — ^ 00), where the non local propag ators G''^^^{i,j;ti,t2,t3,t4) 
represent the corrections at order 0{1/N) to the MF correlations. By an 
appropriated base change, from the unidimensional spin variables Si to the 
eigenvectors r/" of the Jij matrix, we will manage to decouple the interaction. 
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We define r^": 

j 

where a = 1 N and A° is the a-th eigenvalue. The properties of the eigen- 
values and the eigenvectors of the random matrix Jij are the following (see 

i): 

the shape of the eigenvalues density cr(A), for symmetric matrices with random 
and statistically independent elements, such as Jij, in the limit N oo {N 
dimension of the matrix), is: 



a(A) 



2^(4-A2)V2 |A|<2 
|A|>2 



the eigenvectors are statistical variables which have the components indepen- 
dent and following a gaussian distribution defined by the moments: 

^=0, (9) 
(2i^- 1)!! 



ivry^' = ' ' -^ (10) 

they are found to be orthogonal and we can choose them to be normalized: 

N 

Y^rdvH = S^^, (11) 

i=l 

N 

finally, the eigenvectors are uncorrelated among themselves (apart from the 
orthogonality constraint) and they are not correlated to the eigenvalues. 
If we write Sj = y^Tjf in the base of the eigenvectors we obtain: 

= \ E(-/^^" + ^o)(>^")' + f E y'Y^y'Y'Y.vtv^vH- (13) 

We can evaluate the Green functions for the component and, by using (||) 
and (pD-(p!2D, go back to the correlation functions of Si field. 
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We can start evaluating the correlation functions in the static theory. In 



the time-independent limit for the Hamiltonian (13) we may consider the non- 
linear interaction as a perturbation and make a series expansion in the coupling 
constant g. We can show that only the diagrams considered in the Hartree-Fock 
approximation give relevant contributions to the free theory, in the thermody- 
namic limit. This is due to the relations (^-(|1^) satisfied by the eigenvector 
rji. In fact, in the MF limit the relevant contribution to the interaction term 
in the (0) is J^iiv'^Y iv^Y in this way one can see easily that the 

eigenvector index {a for F") plays the same role of the component index in 
the theories of vector fields, where the Hartree-Fock approximation has been 
demonstrated valid when the number of the components goes to infinity. We 
can thus resum all the diagrams at the next orders in g and find a renormal- 



ization of the mass term (the coefficient of the quadratic term in (13)) that for 



r+ is 

2 /T-Ta' /r^2 



and a renormalization of the coupling constant g, that in the same region is 

gr = (2m2)i (15) 

The static susceptibility shows a divergent behaviour for T with the 

critical exponent 7 = 1, in agreement with the results obtained with the replica 
method Q] 



_ 1 v^, ,2 _ - A^) 



^ oci. (16) 



(2m2)V2 r 
The generic 4-point function 

[{SiSkSiSk)^]j = [{{Sif)^{{Skf)^]j + 2[{{SiSk)^f)]j + [{SiSkSiSk)^conn]j (17) 
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is 0{1/N) order and is regular in the critical region. In fact, for T Tc the 
first term is regular and the divergence of the second term 

is compensated by that of the last one 




Also in the replica approach one can verify that the 4-point function (|T7D is 
not singular for T ^ T^. 

In the dynamic case we have to solve the Langevin equation for the time- 
dependent component V^it) (we put Fq = 1): 

= (/5A" - ro)r" - g J2 Y^Y^Y^ ^ vTv^vlv? + C , (20) 

If we define G"(t - 1') = l/2(r"(t)r(f)) and C"(t - t') = (F"(t)F"(t')) we 
find the formal solution 

F"(t) = f' G''{t-t')C{t')dt'+ 
Jo 

-g f dt'G-{t - t') Y: Y^{t')Y^{t')Y\t') Y: r^rirflrfl. (21) 

This is a self-consistent equation, which can be solved with an iterative proce- 
dure. At finite order in g, Y^ is a polynomial in the variables ^"(t). We can 
obtain the dynamic physical quantities averaging G" and C" over the distribu- 
tion of the eigenvalues (H) and of the eigenvectors defined by (i|)-(|ll)- We can 
show that, as in the static case, only the diagrams of the Hartree-fock type are 
relevant in the correction to the Green functions. The correction terms do not 
change the dynamic behaviour of the 2-point functions for T = T^, because 
the time dependent part of the self-energy is regular at T = Tc After some 
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algebraic calculation we find only a mass renormalization effect for C{uj) and 



Gr{uj) 



Cr{uj) 



2 + ATc + r - 2^2{m? - ioj) , (22) 
dXa{X) 



2 ■ 4 



(23) 



^2{m'^ - iuj) + ^2(m2 + iuj) 

where m? is the renormalized static parameter (14), AT^ = — is the 
difference from the renormalized and the bare critical temperature (T° = 2), 
while as usually t = T — T^. The relations (p2|) and (^) are in agreement with 
the critical behaviour indicated in Q by Sompolinsky and Zippelius. The 
susceptibility, according to the definition (^, results 

/■ dX (t(A) 

It is clear that xsci^^i, ^^2) has a finite limit for a; — >• when T > Tc, while at 
T = Tc it shows a critical behaviour such as 

1 , , 

X oc (25) 

For the generic 4-point function the connected terms occur. For example, 
for (^j(ti)sj(t2)sfc(t3)^fc(t4)), we have the diagram of Fig.[l], where, as usual, 
a line with an arrow represents a response function(the time order follows the 
arrow) and one with a cross a correlation function. The diagrams that we have 
to sum in order to evaluate the renormalized coupling constant are drawn in 
Fig. [2] and in the low frequency limit we have: 



= -y2^2(2m2 - iuj). (26) 
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Therefore the total contribution to {C,iSiSkC,k) function, represented in Fig. [3], 
is the sum of the following terms: 



1) --V2 (J2{2m? - 



4 



1 1 



16 



(27) 



2) -iV2(^y^2(2m2 + 



y^2(m2 — icja) + ^2{m'^ + ico'4) 
1 1 
^J2{m? + iuji) + y'2(m2 + ia;2) y'2(m2 - icui) + .y/2(m2 + ztus) 

16 

3) ^ Qv^^2(2m2 - icj)^ Qv^^2(2r772 + 7cu))x 

2i^ I I I \ 4 



2{w? — iuj) ^^2{rn^ + 7u;) / \j2{m? — iuJi) + ^^2{rrfi + zu;2) 

4 

(29) 



/2(m2 - icjg) + ^2(7772 + icU4) 

In the same way we can evaluate all the 4-point functions. The static limit 
derived from the dynamic functions coincides, order by order in with the 
static results. This is true also for the renormalized functions. 

It is not easy to extend this formalism to the short-range model, because in 
that case the eigenvalue density cr(A) is not given by the simple expression (H). 
Although the results of this section have been useful to us for understanding 
some properties of the correlation functions in the critical region, to go beyond 
the MF approximation we have to leave this approach. In the next section we 
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will use the functional integral formalism for a dynamic SG model in finite 
dimensions. The results obtained in this section will be used as comparative 
terms for the long-range limit of the short-range correlation functions. 



4 Fluctuations in short range spin glasses. 
Functional integral formulation. 

To solve a theory with quenched parameters we can use the functional integral 



formulation for dynamics, introduced by De Dominicis |T3[ . So we consider an 
auxiliary field Si{t), [T^, and we define a two-component vector field 
0" = {iSi,Si). One can show that in this formalism the factor isi{t) in a 
correlation function acts like where hi is the external field: 

{smm)Lis,, = ^-^J^^ = Git - (30) 

So isi{t) replaces the noise ^i(t) to generate the response functions. 

After averaging over (in this case we can avoid the replica trick because 
the functional Z is normalized), a 4-spin interaction is generated. It is conve- 
nient to decouple this interaction by a gaussian transformation 0]. Then the 
theory can be defined with a generating functional in the Ql^{t,t') variables: 



2=f n W(l,2)]exp (-/rfld2^Krigf(l,2)A"'3^^gf(l,2) 

+ \nJ[ds][ds]exp{Li{s,s,Qf)'^ , (31) 



where K~j^ = —K^^ [Kij = 1 if i,j are nearest neighbors and zero otherwise) 
and is such that A^'^^'^ = A^^" = A^^^'^ = A^^ai = i and otherwise. 

Li is: 

L, = Lo + J dld2 (12)C(l)^f (2), (32) 
11 



where the field ijjf is defined: 

J I0r = (s,,2s,), (33) 

and Lq is the local part of the theory: 

Lo = JdlJ2 (^-T^'dis,{l) - ros,(l) - + . (34) 

The value of the Jacobi determinant J, associated to the integral formulation 
of the 6 function depends on the discretization chosen to regularize the 



Langevin equation. With the following regularization: 

s{t + 6) - s{t) + e (-^^) = (35) 

where Dl{t) = ft^^ ^i{t') dt', one has J = 1 (and it has been omitted in all the 
previous formula). 

In the previous expressions and in the foUowings of this section we shall 
write 1 for ti, 2 for ^2 and so on. 

The solution of the MF theory is well known [Q, and it is consistent with 
the results obtained in the previous section. Let us consider the fiuctuations 
around the saddle point value Qf^{l,2) = (5"'^(1,2) + 5Q"^(1,2). As a re- 
sult the model can be formulated in terms of the dynamic fiuctuation field 
5gf(l,2): 

Z=f W D{5Qf}exp(L(.,S,Qf + 5Qf)), (36) 

a,/3=l,2 

where L contains the quadratic and the cubic term of the series expansion 
around the MF value: 
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1,2 i,j 

+ In E'^Qf (l'2)C"^"'(l>2,3,4)5Qf(3,4) + 

^ 1,2,3,4 i 

+ ^ E EC"^"'^^(l,2,3,4,5,6)<5Qf(l,2)5Qf(3,4)5Qr(5,6). 

1,2,3,4,5,6 i 

(37) 

The higher order terms can be neglected because we are interested in the 
behaviour of the Green functions near the critical temperature in the param- 
agnetic phase. The C°^t^(1, 2, 3, 4) and C'"^y^^"'{l, 2, 3, 4, 5, 6) vertices are the 
4-spin and 6-spin correlation functions of the one site MF theory described by 
the partition function 

Z, = J[ds,][dsi]exp Lo + y'dld2EQf (l,2)C(l)^f(2) , (38) 
connected with respect to pairs. For example for C"^'''^(l, 2, 3, 4) we have: 

C"^^^(l, 2, 3, 4) = (^ni)^f (2)^7(3)^f (4))Mi.-(C(l)^f (2))A/F(^7(3)^f (4))mf 

(39) 

The form of the functions 0°"^^^ and C°'^^^'^^ is crucial in the following. 
We want to study the universal behaviour of the system near the critical fixed 
point, so we are interesting in the singular part of these functions for T = T^, 
and for u —>■ 0. Thus we consider T = Tc from the beginning in the MF 
theory described by the functional ( |55| ) for the soft-spin model. Because of 
we are not able to compute the 4-point and 6-point functions analytically in 
a closed simple form, a perturbative approach will be used. We perform an 



expansion in the quartic vertex 



using the MF expressions (|22|) 



and ( p3|) for the propagators [{si{t)isi(t'))]j and [(sj(t)sj(t'))] j respectively 

i[{trs^it)ts,{t'))]j = 0). 
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We can demonstrate that the critical behaviour of these functions is deter- 
mined only by the zero loop contributions. In fact the loops that we can form 
with the quartic vertex of Lq, in the uj space, are: 



The first is infrared converging and the latter is only logarithmic infrared 
diverging. In the evaluation of any renormalized correlation function, these 
diverging loops occur always multiplied by correlation functions which are less 
singular (for T T^) than the tree level ones. For example, for the function 

C^^'^^{uJi,U2,uj3,U4) we obtain: 

, UJ2, UJ3, UJ4) = G{uJ4)G{uJi)6{uJ4 + UJ2)S{uJi + ^3) + 

g ■ G{~uj2)G{-uJs) [C{uJi)G{uj^) + G{uj4)G{uji)] S{uji + + cjg + ^^4) + 

■ G{uJi)G{u2)G{u3)G{uj4) ■ \og{u)6{uJi +UJ2 + UJ3 + W4). (40) 



From the expression ( p2D and ( [23| ) for the response and correlation functions, 
we can see that the 1-loop contribution to the connected part of (7^^^^ is 
negligible (it is of order log(co') in the limit uj —>■ 0) with respect to the zero 
loop one (that is of order 1/(0;) ^^^in the same limit). In the same way, we find 
for (7^^^^ that only the connected part is different from zero and is 

gG{-uj2)G{-uj3)G{-uJi)G{uji)6{uJi + U2 + UI3+ '^4)- (41) 

It is reasonable to assume that this phenomenon, which we have seen in a per- 
turbative expansion in g, holds beyond the perturbative theory, so, in order to 
study the critical behaviour of the system we can neglect the loop contribu- 
tions. In these correlation functions were not considered at the critical point 
and a factorized functional expression was proposed for them in the hard-spin 
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limit. The functional form that we obtain is obviously different, but we have 
not yet verified if we obtain also a different value for the universal physical 
quantities (i.e. critical exponents). 

The correlations of the 0° fields, which we are interested in, are related to 
those 5Q'^^ by the relations: 

[(C(l)0f(2))d. = 2Y^{k-%,(af{l,2) + (l,2))^(Q„,)), (42) 

3 

[(C(l)0f (2)0^(3)0^,(4))^], =4 ^(i^'-i),,5:(i^'-i)H 

{ (of (1, 2) + 5Qf (1, 2)) (of (3, 4) + 5Qfi3, 4)) )z.(q.,)+ 

-2{K~^)ikA"'^''^6{l - 3)5(2 - 4) . (43) 

In the next section, therefore, we will to evaluate the correlation functions of 
the fields 6Qf (1,2). 

4.1 The Propagators 

Let us consider the expression ( pTj) with the vanishing cubic interaction. The 
generic propagator 

G-^-^'iz-j; 1,2,3,4) = [(0r(l)0f (2)0j(3)0j(4))^], - [(0r(l)0f (2))d,[(0](3)05(4))^], 
= 4 j:iK-%k{SQf{l, 2)5Qf (3, 4))^(q^^) - 2{K-%A-^^''6{1 - 3)6(2 - 4) , 

I k 

(44) 

is calculated in free theory and will be used to evaluate the corrections of terms 
in the loop expansion. In the free theory G^'^'^'^ii — j; 1, 2, 3, 4) is the solution 
to the following integral system: 

^ (-2 K-\k)A''^''^S{l - 3)5(2 - 4) + C"^^'^(l, 2, 3, 4)) x 

3,4 
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(lQfS^.u^J^. 3^ 4, 5, 6) + ^K-\k)A''^^"'5{3 - 5)5(4 - 6)' 

= -S^,,6p,,K-\k)6{l - 5)5(2 - 6) , (45) 



For each value of /i and u (|45|) is an integral system of four coupled equa- 
tions in the G^^f\ G^^^^^, 0^"^^"", G^^^^ variables. Obviously we consider 
G°^^'^(z, j; ti, ^2, ^3, ^4) invariant under permutation of the index /x, z/ and 
of the relatives times ti, ^2, ^3, ^4- First we solve the system for /x, z/ = 1, 2. We 
are not able to solve the system (^Sf ) exactly, so we use a recursive procedure 
by considering C^^"'^ as a perturbative series in g. For (7 = the coefficients 
C"'^^^(l, 2, 3, 4) are factorized in the time-difference and the integral kernel 
show a complete separation of the internal time. This is the limit in which the 
propagators have been computed in In Fourier space we obtain: 

G]^^\k-Ui,U2,u^,u;^)=Q , (46) 



gr(fe;^i,^2,^3,^4) = K-\k)GMG{-.,)5{u,+.,)5{., + u,) 



Gq {k;uJi,uj2,uj3,uj4) 



G{UJ3)G{—UJ4)X 



K-\k) - GMG{uj2, 



(4J 



K-^{k)-G{uJs)G{-uj,) 

Also at the zero order of the perturbation series in g we consider the mass 
term (m^) renormalized by the interaction. 

At low frequency and for T Tc we have 
K-\k) = -^K-^k) = (4 + AATc + 4r)(l + ck^), G{u}) given by the (0) and 
C{uj) by the (^) and the expression (|^) and (^Sf ) become: 



Gr(^;-.-.-3,^4) = , 4 5(.,+.,)5(0.2+0.3) ^43) 
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6{uJi + UJ3)6{UJ2 + UJ4) + 6{uJi + UJs)6{uJ2 + UJ4) 

4 

X 



X 



(cA;2 + ^2(m2 - io^s/ro) + y'2(m2 + icu^ro)) 
4-2-2 



(50) 



/2(m2 - iujs/To) + p{w? + icus/ro)) 

In the same way, if we calculate the system for i/ = 2,2 we obtain the 
propagator ^2222 ^^isx for = results: 

5{uji + uj^)5{uj2 + 



Gf'^'^{k]UJi,UJ2,0Ji,UJA) = 



4(cA;2 + ^2(m2 - icui/Fo) + ^2(m2 - ic^s/ro)) 
2 



X 



+ 



_(cA;2 + ^2(m2 + ic^g/ro) + ^J^iw? - iui/V^)) 

2 

(cF + ^2(m2 - ia;3/ro) + yj2{m? + ia;4/ro) 
2-4-2 2-4-2 



X 



2(m2 - iu;2/ro) + sj2{w? + ic^s/Fo)) (^2^ - iu^/To) + ^^2(7^2 + ic^s/ro)) 

1 



X 



(cA;2 + ^2(m2 - iu^/To) + ^2(m2 - ia;4/ro)) 

plus another term of the same form proportional to {5{uji + uji)5{uj2 + oj^))- 
For X5.G.(^,<^) = G'^^^^{k;u; — Ui — 102) we obtain the following scaling 



(51) 



behaviour for T ^ T^: 



(52) 



where ^ is the correlation length and rj and z are the usual critical exponents 
that, at the MF approximation, take the value of and 4 respectively. 

At the first order in g we consider C^^-^^ = Cf^^ + Cf^^. For Cf^^ 
we use the former solution and we solve the system in the Gi variables. We 
apply this procedure in an iterative way and we can show that, defining from 
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a diagrammatic point of view 



Gf^'^{i- j]UJl,UJ2,UJ3,UJ4) = Gf^'^{i- j;UJi,UJ2) = 

we obtain, for the following orders in g, recursive expressions that can be re- 
sumed. For istance, the terms that occur for the function G^'^'^^{i—j] uJi, U2, ^3, UJ4) 
can be represented with the diagrams in Fig. [4]. 

The expansion that we obtain can be seen as a perturbative series in the 
quadratic vertex 

G^Z.ih 2, 3, 4)5Qf (1, 2)5Qfi3A). (56) 

In fact, we evaluate the free propagator Gq considering only the disconnected 
part of the C"^^'^'^(l, 2, 3, 4, ), and we use a perturbative approach to take ac- 
count of the connected part. The number of the topologically different vertices 
of form 

CcfZ'(l,2,3,4)<5gf(l,2)5Qf(34) 

are 4 (the number of the superscripts equals one or two, and (7^^^^ = 0). By an 
explicit computation, one can see that the diagrams obtained with the vertex 
( |56| ) are the most diverging ones. 

Indeed, for example, the contribution to the function {6Q^^{1, 2)SQ^'^{3, 4)) 
at the first order in g from the vertex (|56D is 

{6Q'\1, 2)6Q'\3, 4)6Q'\5, 6)6Q'\7, 8)G!!Z{5, 6, 7, 8)) oc 
{SQ'^l, 2)SQ''i5, 6)) {5Q''{3, A)SQ'\7, 8))G^Zi5, 6, 7, 8) 

-^.-^.g.l = g.(l/co') (57) 
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as one can verify from (^) (|23| ) and (^1]), while from the vertex 

C^Z.i5,6,7,S)SQ'\5,6)6Q^\7,8) 
one obtains a weaker singularity of order: 

{SQ^'il, 2)SQ''{3, A)5Q^H5, 6)5Q'\7, 8)C'JZ{5, 6, 7, 8)) cx 
{SQ'\l, 2)SQ'\5, 6)) {5Q'\3, A)SQ'\7, 8))C^Z{5, 6, 7, 8) 

j7.-j7.-9-j7.=9-il/u^'^')- (58) 

At this order in g we do not have any other contribution to {6Q^'^{1, 2)6Q^^{3, 4)), 
because the corrections, that one can obtain with the vertices 

Cnn.(5,6,7,8)5Q2i(5,6)5Q"(7,8) 

and 

ClZ.i^,QJ,8)5Ql\5,Q)6Ql\7,8), 

involved the propagator {6Q^'^6Ql^) that are vanishing at the zero order in g. 

Similarly, other functions can be calculated to verify that stronger singu- 
larity all arise from the vertex 

Therefore, as usual in the case of expansion in a quadratic vertex, we can 
resum the series. Moreover we deal with an order parameter depending on two 
times and we have to consider, in the Fourier space, the integral over the free 
internal frequencies. 

To evaluate the complete propagators at T = T^, we can define the renor- 
malized coupling constant 

gr = g- g'lico, k) + g^l\uj, k) + = , ^ , (59) 

I + ql{uj,k) 
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where I{u!,k), correspondent to the loop 



is given by the integral: 
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Ii{uj,k)= I — - — X 

1 16 



(cfc2 + ^2iuj/VQ + ^-2i{uj - uj)/Vq ^-2iuj/VQ + ^2iuj/Vi 



:X 



' { TT-^l • (60) 



cA;2 + ^-2iuj/To \(^cu/ro) 



/ cfc2 \ 

Fi is a function of the variable _ , that exhibits a constant hmit for 
^ and for ^ . The couphng constant in the low frequency and small 
moments limit is, according to (EUf) and to 



g, = {ck' + ^-2iQlT^)l/F^. (61) 

To evaluate the total contribution of the loop corrections we must consider the 
term given by the following diagram: 



with renormalized vertices. The value of the loop is: 

duj 1 



k) 



1 



27r cA;2 + y/-2zcj/ro + ^ -2i{u) - uj)/To ck^ + ^2iuo/To + ^2i{u: - uo)/Vq 
1 1 



+ 



cA;2 + J-2iuj/VQ + J2i{uj - uj)/Vq ck^ + j2iuj/To + J-2i{iu - uj)/T, 



X 



2iiu/To + \/2iuj/To \/-2i{uj - iu)/To + ^2t{iu - iu)/To 

„ / cA;2 \ 



X 



where -^2(2;), like -Fi(x), has a constant limit for — > and k ^ 0. 
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The connected part of the free propagator G^^^^ aX T = is computed 
adding up all the diagrams in Fig. [4]. We have: 



k^ + 2iuj/To \ck^ + J2iij/To 



+Gl^^\k;uJi,uj2)Gl^^\k;uj3,uJi)gr{k;-ij)^ + 

+Gl^^^{k; uji, LJ2)Gf^^{k] Us, ujA)gr{k] uj)gr{k; -uj)-^ 

^1 

6{uj1 + uj2 + uj3 + ujA) 
(63) 

which corresponds to the diagrams in Fig. [5] with the constant gr given by 
( |HT| ) and lD = LJ2 — i^i- 

5 Conclusion 

The connected term of the propagators is zero only in the limit of two complete 
time separations. In all the other cases we must calculate the complete corre- 
lation function. In this way we compute G^^'^^ for all the values of the indices 
a(3'^5 and for all the time distances. It is easy to see that the long-range limit 
{k 0) of the connected part of the expressions (|63|) coincides with the sum 
of the terms (^7]) (pSf ) ( ^Of ) of the third section. On the level of the Gaussian 
approximation, i.e. cubic interactions are neglected, the connected part of the 
propagators does not contribute to the susceptibility and the dynamic scaling 



([52|) is correct. 

Therefore, we have analysed the critical behaviour of the propagators of the 
soft-spin model in the quadratic approximation and we have put the bases for 
a short-range theory of SG in the renormalization group formalism. In fact the 
expressions that we have derived could be used to evaluate the contributions of 
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the Feynmann diagrams that occur in the loop expansion when the cubic term 
of the Lagrangian (|3^) is considered non vanishing. The expression which we 
have obtained in the soft-spin case are quite different from those obtained by 
Zippehus in ref . . It is certainly interesting to understand if the value of the 
dynamical critical exponent is affected by this difference. The computation of 
the loops will be crucial to clarify this point. 
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CAPTIONS FOR ILLUSTRATIONS 



Fig. 1: Diagrams for the correlation function [{^iSiSk^k)] 
in the series expansion in g. 

Fig. 2: 1-particle irreducible diagrams that contribute to 
the renormalization of the coupling constant. 

Fig. 3: The renormalized function [{^iSiSkik)]conn- 

Fig. 4: Diagrammatic representation of the propagator 

^1221^-^ — j; a;2, a;3, a;4) in the series expansion in g. 

Fig. 5: The connected part of the propagator 



26 



